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Abstract 
Introduction. Some of the major factors in the occurrence of vibration of units, causing fatigue failure of the housings 


of elements, pipelines, and failure of pump elements, are pulsations of the working medium in the hydraulic systems of 
machine tools, fuel feed systems of aircraft engines and liquid-propellant engine supply. This study aimed at the 
implementation of a method for determining the dynamic characteristics of a volumetric pump using special bench 
systems, and the comparison of the calculation results to the experimental data. The stages of calculating the dynamic 
characteristics of a volumetric pump were described, taking into account the pre-developed special bench systems on 
the example of an external gear pump with a capacity of 14 cm*/rev. The implementation of V. P. Shorin's load 
variation method using special bench systems developed by the authors with predetermined dynamic characteristics was 
shown. The main stages of the methodology for determining the dynamic characteristics of a gear pump were described. 
Materials and Methods. Methods of spectral analysis of pulsating pressure were applied in the work. Pulsations 
of fluid flow at the pump outlet were determined using the impedance method, the method of load variation, and 
special bench systems. 

Results. The paper implemented a technique for determining the dynamic characteristics of a gear pump in the drive 
shaft speed range of 500-2500 rpm for four harmonic components of the vibration spectrum in a wide range of dynamic 
loads (from inertial to capacitive nature). The bench systems yielding the calculation of the dynamic characteristics of 
the pump with a minimum error based on the condition of matching the dynamic load and the source of vibrations were 
analyzed. The developed approach to the evaluation of the dynamic characteristics of the pump was verified through 
comparing the calculated and experimental data of pressure pulsations in the bench systems with choke, cavity and an 
extended pipeline at the pump outlet. 

Discussion and Conclusions. The method for determining the dynamic characteristics of a volumetric pump was 
implemented using special bench systems developed by the authors. The research results show that the gear pump under 
study can be considered as an independent source of flow fluctuations, for which the deviation of its own dynamic 


characteristics from the average values does not exceed 10% for the first harmonic component. 
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Introduction. The main disadvantages of pumps are the acoustic noise that accompanies their operation, the 
vibration of the body, connected pipelines, and the resulting structural damage. 

The traditional way of improving pump designs, used in mechanical engineering and consisting in smoothing out the 
uneven flow of the liquid, is hampered by the lack of means of direct measurement of fluid flow at the pump outlet. 
This limits any experimental evaluation of the design efficiency and finishing work on pumps. In addition, this limits 
the analytical description of the pump as a source of vibrations required by hydraulic system developers to provide the 
target dynamic quality of the systems by themselves. 

For this reason, the research on the modeling of dynamic processes in hydromechanical systems [1,2] and the 
application of new approaches to solving problems of hydraulic system dynamics [3, 4] based on physical principles are 
pressing issues. 

Materials and Methods. The implementation of V. P. Shorin's method of load variation [4] is presented. It consists 
in determining pressure pulsations behind a pump interacting with several bench systems, whose dynamic 
characteristics are calculated in advance. Further, the processed frequency-response functions of the pressure using the 
impedance of the bench systems are recalculated into the pump's own dynamic characteristics — pulsation performance 
and impedance. The development of the approach provides that special bench systems are proposed. They take into 
account the connecting fittings, adapters and internal channels of the pump, which makes it possible to more accurately 
determine the dynamic characteristics of pumps, as well as to verify the applicability of the pump model used in the 
form of an equivalent source of vibrations. 

Research Results 

Method for determining the dynamic characteristics of a gear pump 

Samara University has developed a computational and experimental method for determining the dynamic 
characteristics of a gear pump (pulsation performance and impedance), which is based on the research of L. Thévenin, 
E. Norton and V. P. Shorin [5, 6]. The developed methodology includes the following main stages: creation of models 
and calculation of dynamic characteristics of special bench systems; experimental determination of pressure pulsations 
at the pump outlet; calculation of the spectrum of excited vibrations in various bench systems; determination of 
approximating dependences for individual harmonic components; calculation of the pump own dynamic characteristics 
(pulsation performance and impedance); evaluation of the independence (stability) of the dynamic characteristics of the 
pump from the characteristics of bench systems; verification of the applicability of the pump model used as an 
equivalent source of flow (or pressure) fluctuations. The following is an example of the implementation of the proposed 
technique using an external gear pump with a specific capacity of 14 cm*/rev. 

Models and calculation of dynamic characteristics of special bench systems 

Basic options of bench systems at the liquid outflow during pump tests were considered. In concentrated parameters, 
the bench system can implement active flow-dependent, inertial, and elastic loads. In distributed parameters, the bench 
system can be provided in the form of an extended cylindrical pipeline implementing an active frequency-independent 
non-reflective load. 

Pump models in the form of equivalent source of flow fluctuations (ESFF) or pressure (ESPF) are determined using 


constant flow source Q,, an ideal source of pressure fluctuations p (or flow q ) and internal impedance Z, [4]. Figure 1 


shows pump models in the form of an equivalent source of pressure fluctuations (and flow) with an attached system. 


bench 
pump | system pump | system 


a) b) 


Fig. 1. Pump models in the form of an equivalent source of: a — pressure fluctuations; b — flow with attached system 


Q,,P, — average flow and pressure; p,,,g,, — pump's own pulsating performance in terms of pressure and flow, 


Z.,,»Zem.c. —~ impedance of pump and bench system (adapted from [4]) 


ur 
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In this connection, the task of selecting a bench system is to provide such system parameters at which the best 
conditions for high measurement accuracy are realized. The task of selecting the parameters of the attached system 
should be carried out when the impedance changes in the range from -00 to +00, 

Figure 2 shows schematic diagrams of the bench systems. 
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Fig. 2. Schematic diagrams of the bench systems for determining own dynamic characteristics of pumps: a — with active throttle 
load; b — with active non-reflective load; c — with capacitive load: 1 — test pump; 2 — load choke; 3 — flow meter; 
4 — pressure sensor; 5, 6 — connecting fittings; 7 — main line; 8 — extended pipeline; 9 — cavity 


Computational models of impedances of special bench systems: 
— with active throttle load (throttle): 


ons, LE (1) 
. 1 Jet anes, 
Zo (Q,, f) = |(Z,, (Q, ) +ReZ,,, (Q, yy [ann dl e  Zin(Qo)+ReZnp( Qo) 


mp 


— with active distributed load (extended pipeline): 


pa (2) 


— with capacitive load (cavity): 


‘ (3) 
2 | (2af,) LC-1 
Zone (Qos fi) = q|( REZ ses (Qo) + REZyys (Qn) + REZ (Op) + Zip (Qa)) a ; 
k 
wr (27h; y Lc-1 
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The formulas indicate: Z,, (Q,) — impedance of the throttle; ReZ,,(Qo) — real part of the impedance of the 
pipeline, including the connecting line, fittings, internal channels in the pump and _ throttle; 
ReZ nam (Qo), REZryp6(Qo), ReZu.c.(Qo) — real part of the impedance of the fittings and the internal channel in the 
pump and cavity with laminar, turbulent fluid flow behind the pump and local resistances, respectively; f, — 
frequency k — th of the harmonic component of the vibration spectrum behind the pump;  — density of the working 
environment; /, — pipeline length, including the length of main line 6, connecting fittings 5 and 7 and internal 


channels of the pump and throttle (not shown in Fig. 7); S,,, — cross-sectional area of the pipeline; L — “inertia” of 
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the connecting fittings and the internal channel in the pump and cavity; C — elasticity of the cavity; k — number of 


the harmonic component of the vibration spectrum; j — imaginary unit (j = |. 


Using a bench system with a throttle at the pump outlet, three inertially active dynamic loads were implemented, in 
which the active resistance was changed through changing the overlap area of throttle Fo: 

— for bench system no. 1, overlap area Fo, =9.290x 10° m’, 

— for bench system no. 2, overlap area Fy:=7.299x10° m’, 

— for bench system no. 3, overlap area Fo3=2.247x10° m’. 


Using a bench system with a cavity at the outlet of the pump with a volume of 2 liters, an inertial elastic dynamic 


load was realized. 


An active frequency-independent load was implemented using a bench system with an extended pipeline 


length of 106 m. 


The dependences of the impedance (module and phase) on the frequency were calculated according to the models of 
bench systems using formulas (1)—(3). Table 1 for the Ist harmonic shows the regression functions of the impedance 
module of the bench system Z;(f) and phases B;(f) for the analyzed bench systems. 


Table | 


Regression functions of dynamic characteristics of bench systems according to the 1st harmonic 


No. of stand Bench system impedance 
De Hensnaytomee module Bench system impedance phase §; (f), 
according to Z,(f) ke rad 
Figure 2 mas? 
a—no.1 ainciioat 1.412-10’f + 2.201-10" 7.066: 10™f + 0.199 
a—no.2 a 1.610°10'f + 4.25810" 6.906- 10°F +0.09039 
a—no.3 4.980-10'f + 1.503-10” 2.723-10°f + 0.251 
.308-10"f? — 1.416-10’ 4.234-10° f +1.495-10° 
caw With cavity at pump 3.308:10°f — 1 ale 10°'f+ Grete - f ; 
outlet 1.873-10 1.430-10° f +2.930-10 
With extended 
c —no.5 pipeline at pump 1.897-10'° 0 
outlet 


Experimental determination of pressure pulsation at the pump outlet 


The background information for determining the dynamic characteristics of the pump was the measurement of 
pressure at the pump outlet. An example of an oscillogram and a pressure pulsation spectrum at the pump outlet in the 
bench system with an active distributed load is shown in Figure 3. Operating mode of the pump in this case is as 
follows: average pressure PR, = 23.8 MPa, drive shaft speed n=1000rpm. 
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Amplitude of pressure pulsations, MPa 


167 334 501 668 835 1002 1169 1336 1503 1670 1837 
Frequency of oscillation spectrum components, Hz 


b) 
Fig. 3. Pressure at the pump outlet: a — oscillogram; b — spectrogram 


The amplitude-frequency response of pressure pulsations at the pump outlet in the speed range of the pump drive 
shaft n=500—1250 rpm is shown in Figure 4. The regression functions built in the Microsoft Excel environment using 
the Trendline tool are plotted with a solid line. 
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Amplitude of pressure pulsations, MPa 
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Frequency of vibration spectrum components, Hz 
———— 3 
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i: -_ 
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Fig. 4. Amplitude-frequency response of pressure pulsations at the pump outlet by individual harmonic components 
(dots indicate data on: ® — first, ® — second, ® — third, @ — fourth harmonics, lines indicate regression functions) 


The analysis of the oscillogram, spectrogram and amplitude-frequency response of pressure pulsations showed the 
following: the type of process was polyharmonic steady; the number of recorded harmonics was 8, analyzed! — 4; the 
amplitude range of pressure fluctuations was 0.28—0.52 MPa; the drive shaft frequency range was 500-1250 rpm, the 
analyzed harmonic components of the spectrum 83-833 Hz. 

Fourier series expansion was used for harmonic analysis of periodic signals [7]. To decompose into a Fourier series, 
sections of the pressure oscillogram measured relative to the reference signal at the rotor frequency were used, i.e., one 
section of the oscillogram was extracted for a full rotation of the rotor. The reference signal was recorded at a fixed gear 
position, set at the beginning of recording the pressure waveform for five bench systems. In this way, the amplitude- 
frequency (FRF) A;,(w) and phase-frequency (PFC)@;,(w) characteristics were determined. 


' Tt is assumed in the work that the extreme highest analyzed harmonic in amplitude should not be more than 5 % of the amplitude of the first 
harmonic. 
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Calculation of the spectrum of excited vibrations in various bench systems and determination of approximating for 
individual harmonic components of pressure fluctuations 

Based on the experimental determination of pressure in five bench systems: with a throttle (no.1—3), with a cavity 
(no.4) and an extended pipeline (no.5) at the pump outlet, the spectra of excited vibrations were calculated. Table 2 
shows the regression functions of the amplitudes A,(f) and the initial phases of pressure pulsations for the 1st harmonic. 


Table 2 
Regression functions of amplitudes and initial phases of pressure fluctuations 
No. of stand sas 
system Pressure pulsation amplitude Initial phase ph Frequency 
. Bench system type kg pressure pulsations 
according to Api &), | tad range 
Figure 2 be 
—no.1 1.003-10°f 0.543 
a With throttle at : 
a—no.2 1.243-10°f —0.280 
pump outlet x 
a—no.3 6,129-10°f —5.941 
With cavity at pump 1.343-10°F — 4.298-107f + 83-417 
poe outlet 3.913-107 a 
With extended 
c —no.5 pipeline at pump 2.694: 10°F —0.658 
outlet 


Calculation of the dynamic characteristics of the pump 
To calculate the dynamic characteristics of the pump (impedance and pulsation performance of the source), it is 
required to conduct an experiment with at least two dynamic loads”. Therein, according to the model of an equivalent 


source of flow fluctuations refined by the authors using analytical dependences p, ( f ) and Z, ( f ) from Tables 1 and 2, 


the pump impedance is calculated from the formula: 


jon 
P, 
L =) 4 (4) 
PL 
P, Z, 
where Z, — pump impedance; ( Z, = [Zee ) — pump impedance module; (|Z,|=./ReZ,,° + ImZ,," ), ReZ, ,lmZ,, — 
real and imaginary part of pump impedance; @, — pump impedance argument; Z,,Z, — input impedance of dynamic 


loads no. 1-2; (Z, =|Z,|e", Z, =|Z, |e"), 2,8, — impedance arguments of dynamic loads no. 1-2; p,,p, — 


pressure pulsations behind the pump in bench system no. 1-2; ( p, = A,,e"", p, = A,,e’” ), A,,,A,. — amplitudes of 


flow pulsations behind the pump in bench system no. 1-2; g,,g, — initial phases of pressure pulsations of dynamic 
loads no. 1-2. 

Then, using the analytical dependences from Tables 1 and 2, the pulsation capacity of the pump was determined 
according to the model of an equivalent source of vibrations refined by the authors (variable component of volumetric 
flow q, and pressure p, ): 


ma 14 
P2 Z, Z, 
uo ? Py = (5) 
Z, l P2 PZ, 
P P, 2, 


? For a sufficient sample size, it was required to use at least four dynamic loads Z; , that form six non-repeating (including i # j) combinations: Z,, , 


Zi3» 2g Za3» Loy, Lrg 
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where q,,p,, — variable component of the volume flow (gq, = Ae") and pressure ( p, =A,,e""") of the pump; 


A,,»A,,— pump flow and pressure pulsation amplitude; g,,~, — argument of the variable component of the volume 
flow and pressure of the pump. 
Of the calculated characteristics according to ESFF (q, (f)) and ESPF ( p,(f)) models, using formula (5), only 


one will be stable, since the equivalent oscillation source has an independent dynamic characteristic only for one of the 
parameters (pulsating flow or pressure). 

The calculation of impedance Z, and the variable component of the volume flow g, of the pump was carried out 
according to formulas (4)—(5) for five dynamic loads, i.e., using ten different combinations of these loads: “1_2”, “1 3”, 
“237,14, 24, 34 “15”, 25”, “35”, “45” (Table 1). 


For further calculation and analysis, a characteristic of the pump flow pulsation amplitudes A, was obtained from 


ten flow pulsation amplitude curves. The values of the amplitudes A, (at each frequency) were calculated using the 
formula from source [8]: 

Aw = Au £ BA as 

where Au — average amplitude of pulsating flow; AA,, — confidence interval of the amplitude of pump flow 


pulsations (AA, = aes, o — standard deviation of the amplitude of the pulsating flow rate from the average value; 


qu in 


m — number of dynamic load combinations (m=10); t— Student's coefficient (t = 2.262, was calculated for the 
confidence level 0.95). 


The confidence interval of amplitudes AA, was formed from errors (instrument and method) and the degree of 


dependence of the dynamic characteristics of the pump A, ( f ) on the bench system. The instrument error consisted of 


the error of determining the pressure and phases, the method error included the error of the unaccounted influence of the 
dynamic properties of the drive mechanical system. 

The phases of flow pulsations gy, and their deviations were calculated in a similar way. When calculating the 
phases, combinations “1 2” and “1 5” were excluded due to the degeneration of the variable component of the volume 
flow of the pump g, when using close values of characteristics at the same frequency: 9, ~ @, and f, = f, . 


The results of calculating the intrinsic dynamic characteristic of the pump in the form of pulsation performance and 


impedance are shown in Figure 5. Phases, for convenience, are presented in degrees. 
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90 
s 
= 60 — 
2 
eS 30 
g 
= oO 
Pes) 0 
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gE 
2 a -30 
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Ss 
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Frequency of harmonic components of the vibration spectrum, Hz 
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Pump impedance phase, deg 
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83 264 445 626 807 988 1169 1350 = 1531 


Frequency of harmonic components of the vibration spectrum, Hz 


d) 
Fig. 5. Dynamic characteristics of the pump: a — flow pulsation amplitude; b — flow pulsation phase; c — impedance modulus; 
d — impedance phase (solid lines indicate data on: eum — first, @======—" — second, 
eum — third, emmmmme — forth harmonics) 
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The amplitude of the flow pulsations at the 1st harmonic has a monotonically increasing character. With an increase 
in the harmonic number, the slope of the flow pulsation amplitude curve decreases to an almost constant value to the 4th 


harmonic. The values of amplitudes A, at the higher harmonics are by an order of magnitude less than at the Ist 


harmonic. 
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To evaluate the realized dynamic loads, the criterion of matching the impedance module \Z,,| of the oscillation 


was used [99, 10]. The impedance modules should 


cm.c 


source and the input impedance module of the dynamic load |Z 


be comparable (|Z, ), so that the quantitative values of the dynamic characteristics of the pump were 


=|Z 


cm.c 


u 


u 


was used. At 


determined with minimal error. For comparison, ratio ~1, the dynamic load was considered 


cm.c. cm.c. 


consistent with the source. 
Due to the difference in the curve shapes of the pump impedance modules and dynamic loads in the frequency range 


u 


of 83-417 Hz, it was rational to use average ratio , that provided considering the contribution of most of 


cm.c. 
min **Fmax 


u 


the frequency range, i.e., ratio according to the Ist harmonic would be: 


cmc. 


Finin -+-Simax 

— for active throttle load no.1 — 1.16; 

— for active throttle load no.2 — 0.8; 

— for active throttle load no.3 — 2.65; 

— for capacitive load no.4 — 4.64; 

— for active distributed load no.5 — 0.35. 

That is, the dynamic loads of bench systems no. 1, 2 were the most consistent with the source of vibrations. Thus, 
the dynamic characteristic of the pump, calculated using dynamic loads of bench systems no.1, 2 (“1_3”, “2.3”, “1_4”, 
“2 4”, “1_ 5”, “2 5”) was determined with a minimum error. 

Evaluation of the stability of the dynamic characteristics of the pump from the characteristics of bench systems 


Stability of the dynamic characteristics A, and 9, (IZ, and @, ) of the pump was evaluated using mathematical 


Statistics tools [8]. Thus, the stability of the characteristic A,, was estimated through calculating the ratio of the 


confidence interval |AA at this frequency* The 


qu 


at frequency f to the average value of the characteristic [Anat 


resulting relative deviation A, should not exceed the specified deviation A,,,, : 


ae MA, on 
q | | — dong? 
quMe 
where [A uae — average median value of the amplitude of pump flow pulsations; A,,,, — specified relative deviation 


of the pump flow rate pulsation amplitude from the average value. 


With relative deviation A, 0 (when |A4,, —> 0), the stability of characteristic A,, was maximum. Based on the 


experience of other researchers, it is preferable to select value A,,,, within 5-30 % [11-17]. 


For a single-digit evaluation of the dynamic characteristics of the pump, it is rational to use the value of the relative 
deviation. The calculation of such a value is possible through estimating the average median value of relative deviations 


AA 
in the frequency range of each harmonic ( fyi,---fmue): Mgve =—— . Value A,,,, is efficiently used, if 
PN fin Sma 
MA, . . : . . 
< 0.3, i.e., when the confidence interval AA,,, does not exceed 30% of the average amplitude of pulsating 
A 
Fin Sina 


qu * 
The results of assessing the stability of the dynamic characteristics of the pump from the characteristics of the bench 
systems’ were: according to the Ist harmonic — 2 %, according to the 2nd — 8 %, according to the 3rd — 15 %, 
according to the 4th — 48 %. The results of the calculation of the dynamic characteristics of the pump showed that for 


> Stability of other intrinsic characteristics of the pump is evaluated similarly. 
“Tt is required to consider the method error of the approach and the instrument error of determining the pressure. 
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four harmonics in the frequency range of 500—2500 rpm, the gear pump under study should be considered according to 
the model of an equivalent source of flow fluctuations. 

Checking the applicability of the pump model used as an equivalent source of flow fluctuations 

To test the developed approach for evaluating the dynamic characteristics of the pump, a calculated determination of 
pressure pulsations in bench systems with a throttle, a cavity and an extended pipeline at the pump outlet was 
performed. Using the found dynamic characteristics of the pump, the amplitudes of pressure pulsations in the analyzed 
bench systems were calculated, which were compared to the experimental ones. The results of the calculation of 


amplitudes A, ( f ) in a bench system with an extended pipeline are presented below. The amplitudes of pressure 


pulsations p, were calculated in accordance with the ESFF model from the formula: 


qi . Z, . Z; (6) 
| aaa 
Z,+Z, 
The results of calculating the amplitudes of pressure pulsations behind the pump are shown in Figure 6. The dotted 
line indicates the amplitudes obtained experimentally, and the solid line — calculated from formula (6). 


= 
f=) 


0.8 


0.6 


0.4 - 


0.2 


Pressure fluctuation amplitude, MPa 


0.0 


83 283 483 683 883 1083 1283 1483 1683 
Frequency of harmonic components of the vibration spectrum, Hz 
Fig. 6. Pressure frequency response for the bench system with extended pipeline length /=106 m (solid lines indicate calculated data 


On: Gomme — first, Ge — second, eum — third, em — forth harmonics; 


dotted lines indicate experimental data) 


Convergence of the results of the verification calculation was evaluated according to the formula of deviation of the 
difference between the experimental and calculated values from the experimental data. The results of the test calculation 
for the bench systems with a throttle, a cavity and an extended pipeline at the pump outlet showed that according to the 
1st harmonic, the relative deviation from the experimental data for the considered bench systems did not exceed 10 %, 
the minimum relative deviation was characteristic of a bench system with an extended pipeline at the pump outlet — 
2%. The amplitudes of the pressure pulsations of the 2nd—4th harmonics were within the corridor calculated from the 


formula A, =A,+A, (where A, — the average amplitude of the pulsating pressure). 


Discussion and Conclusions. Methodology for calculating dynamic characteristics of the pump in the form of an 
equivalent source of fluid flow fluctuations was implemented on the example of a gear pump with a capacity of 14 
cm’/rev using special bench systems developed by the authors and an updated model of impedance and pulsation flow 
rate. It was shown that the investigated gear pump of external gearing was rationally considered according to the model 
of an equivalent source of flow fluctuations. The proposed approach was verified in the form of calculation of pressure 
pulsation amplitudes in the bench systems and comparison to experimental values. The deviation of the calculated 
amplitudes of pressure pulsations from the experimental values did not exceed 10%. It is worth noting that the impact 
on the dynamic characteristics of the gear pump of the drive mechanical system and the wear of the gear pump parts 
remained outside the scope of the review. Based on the work carried out, it was possible to formulate directions for the 
development of the proposed approach: 

— assessment of the dynamic characteristics of the gear pump during the entire service life; 
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— development of a four-pole model of a gear pump as an equivalent source of vibrations, taking into account the 
mechanical drive and hydraulic output systems; 

— development of a gear pump model as an equivalent source of vibrations, taking into account its design features, 
the displacement process, reverse water hammer in the chamber between teeth, the work of the locked volume, and 
leakage from the injection area to the suction area. 

Here, it is required to expand the nomenclature of the tested volumetric pumps for the entire class of volumetric 
displacement machines to confirm the operability of the proposed approach. 
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